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Motivated by the colossal negative thermal expansion recently found in BiNiCL, the valence transition ac¬ 
companied by the charge transfer between the Bi and Ni sites is theoretically studied. We introduce an effective 
model for Bi-6s and Ni-3d orbitals with taking into account the valence skipping of Bi cations, and investi¬ 
gate the ground-state and finite-temperature phase diagrams within the mean-field approximation. We find that 
the valence transition is caused by commensurate locking of the electron filling in each orbital associated with 
charge and magnetic orderings, and the critical temperature and the nature of the transitions are strongly affected 
by the relative energy between the Bi and Ni levels and the effective electron-electron interaction in the Bi sites. 
The obtained phase diagram well explains the temperature- and pressure-driven valence transitions in BiNiC >3 
and the systematic variation of valence states for a series of Bi and Pb perovskite oxides. 

PACS numbers: 71.10.Fd , 71.30.+h , 75.25.Dk, 75.30.Kz 


Perovskite transition metal (TM) oxides (general formula: 
AB O 3 ) have been providing central issues of phase transi¬ 
tions and strong electron correlations in condensed matter 
physics mi. They exhibit a wide range of novel mag¬ 
netic, dielectric, and transport properties: for example, the 
large negative magnetoresistance in Lai_ x Sr T MnO ,3 i§3, 
the spin-state transition in Lai Sr.,Co 0.3 | 6 . 7j, the metal- 
to-insulator transition in .RNi0 3 (R: rare earth element) (3, 
and the ferroelectric to quantum paraelectric transition in 
Bai_ x Sr x Ti0 3 @,[13]. In these phenomena, the central play¬ 
ers are the electrons in 3d orbitals of the /7-site TMs hy¬ 
bridized with oxygen 2 p orbitals. The A-site cations, on 
the other hand, are usually inert and have been regarded as 
“stagehands”: they control the electron filling and bandwidth 
through their valence state and ionic radius, respectively. 

Peculiar exceptions to the above standards have recently 
been found in several perovskite TM oxides, in which the A- 
site cations play an active role as “valence skipper”. In these 
compounds, not only the B -site 3d electrons but also the va¬ 
lence s electrons in the .4-site cations significantly contribute 
to the electronic properties. In the valence skippers, the outer¬ 
most s orbital prefers closed-shell configurations s° or s 2 , and 
tends to skip the intermediate valence s 1 . This is attributed to 
the effective attractive interaction between s electrons E3- 
131I . and hence the .4-site valence state can be actively con¬ 
trolled through electronic degrees of freedom. Owing to the 
multiple electronic instabilities in both A- and B -site cations, 
the TM oxides with the A-site valence skipper have a potential 
of new electronic phases and functions. 

The colossal negative thermal expansion (CNTE) material 
BiNi0 3 |l4] is one of such candidates; both Bi-6s and Ni- 
3d electrons are expected to play a key role in the large vol¬ 
ume change jl5i 11611 . At ambient pressure, BiNi0 3 has a 
unique valence state, where the average valence of Bi is 4+ 
but it is disproportionated to 3+ and 5+, while the valence 
of Ni is 2+ 01711 : the A-site Bi cation exhibits the valence 
skipping nature. Bi 3+ and Bi 5+ are spatially ordered in a 


checkerboard-like pattern, which can be regarded as a charge 
ordering (CO) at the Bi sites. The electronic and lattice struc¬ 
tures are drastically changed by applying pressure; the system 
exhibits an insulator-to-metal transition with the structural 
change from triclinic to orthorhombic, where the electrons are 
transferred from Ni to Bi and the valence state changes from 
Bi 3+ 0 . 5 Bi 5+ o. 5 Ni 2+ to Bi 3+ Ni 3+ . At the same time, the unit 
cell volume largely shrinks about by 5 % JB Ellll. Un¬ 
der pressure, this phase transition is also observed by raising 
temperature (T), which is termed as the CNTE. 

The large volume shrinkage in the CNTE is attributable to 
the valence increase at the B sites which determine the lat¬ 
tice parameters in perovskites in general. However, the mech¬ 
anism of the valence transition behind it, involving both Bi 
and Ni sites and the charge transfer between them, remains 
to be clarified. This is owing to the fact that, as mentioned 
above, most of the previous studies for perovskite TM oxides 
have been generally focusing on the TM 3 d and oxygen 2 p 
electrons, which usually govern the electronic properties, and 
hardly incorporate the valence state in the ,4-site cations ex¬ 
plicitly. 

In this Letter, we present a microscopic theory for the va¬ 
lence transition behind the CNTE in BiNi0 3 , which treats the 
active electronic degrees of freedom in both A- and /j- site 
cations on equal footing. We introduce a simple but realistic 
effective model for BiNi0 3 , and clarify the ground-state and 
finite-'/' phase diagrams within the mean-field (MF) approxi¬ 
mation. We will show that the valence transition is controlled 
by the relative energy between Bi and Ni levels as well as the 
electron correlation at the Bi sites. The charge transfer be¬ 
tween Bi and Ni is caused by spontaneous symmetry breaking 
associated with a bipolaronic CO at Bi and antiferromagnetic 
(AFM) order at Ni, both of which favor a commensurate fill¬ 
ing in each band. Our result will provide a useful guide for 
further exploration of larger CNTE in related materials. 

Let us construct an effective model for BiNi0 3 . Accord¬ 
ing to first principles band calculations, the Bi-6s and Ni- 
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FIG. 1: (Color online) (a) Schematic electron configurations for Bi 
6s and Ni 3 d e g orbitals in BiNiCL and (b) those in the effective 
model in Eq. Q] CO, PM, AFM, and CT denote charge order, para¬ 
magnetic, antiferromagnetic states, and charge transfer, respectively. 
Schematic pictures of (c) pervskite structure for BiNiOs and (d) the 
projected two-dimensional lattice structure for the effective model. 
The shaded area represents the mean-field unit cell. 

3 d e g bands are located near the Fermi energy, and the O- 
2 p bands strongly hybridize with them 0. Here, we take 
into account the Bi-6s and Ni-3d e g orbitals, while assum¬ 
ing that the role of 0-2 p is effectively incorporated in the en¬ 
ergy levels and the other parameters. The electronic config¬ 
urations in the insulating Bi 3 + o. 5 Bi 5 + o.sNi 2 + 0 2_ 3 and the 
metallic Bi 3+ Ni 3+ 0 2 ~3 phases are schematically illustrated 
in Fig. Q2 a ). For further simplicity, we omit the twofold de¬ 
generacy of e g orbitals as shown in Fig. |TJb); the essential 
physics related to the valence transition will be retained as 
described below. Considering the uniform electron configu¬ 
ration along the z axis in the real compound 0, we adopt 
a two-dimensional lattice of Bi and Ni obtained by projecting 
the perovskite structure in Fig. |TJc) on the xy plane, as shown 
in Fig. Q]d). The Hamiltonian is given by 

Ni-Ni Bi-Bi 

n = t N J2 ( a L a ^ + H - c -) +tB J2 ( b L b vr+ h.c.) 

(ij)cr (ij)cr 

Bi-Ni 

+ fBN ^ ( a L^'cr + H.C.^j 

(+/V 

Ni Ni Bi 

+ A J2 n ^ + UN J2 n it n 4- + Ub J2 

icy icy icy 

Bi—Bi Bi—Ni 

+ V B ^2 nfnf+V B K ^2 n f n f^ (!) 

(ij) {ij) 

where and b lrJ represent the annihilation operators of elec¬ 
tron with the spin <r(=t, 4-) at the Ni and Bi sites of i-th unit 
cell, respectively; nf a = a\ a a i( j and nf a = b} a b i<7 . In Eq. (Q}. 



FIG. 2: (Color online) (a) Ground state phase diagram. Solid and 
broken lines denote the second and first order transitions, respec¬ 
tively. Schematic pictures of charge and spin configurations: (b) 
uniform-I, (c) CO+AFM, (d) uniform-II, (e) sub-I, and (f) sub-II. The 
circles and arrows represent the charge densities and spin moments 
at each site, respectively. 


the first and second lines represent the electron hopping on 
the Ni-Ni, Bi-Bi, and Bi-Ni bonds. In the third line, the first 
term is the energy difference between the Bi-6s and Ni-3d lev¬ 
els, the second term is the on-site Coulomb interactions on the 
Ni sites, and the third term is the effective interaction that de¬ 
scribes the valence skipping nature of Bi iTHT - fpll : we consider 
not only positive but also negative values for U B ■ The fourth 
line represents the intersite Coulomb interactions on the Bi-Bi 
and Bi-Ni bonds. 

We obtain the phase diagram for the model in Eq. £[]) by the 
MF approximation with decoupling the two-body interaction 
terms in the third and fourth lines in Eq. (|T|) as nn' ~ n{n') + 
(n)n' — (n) (n 1 ) . We take the unit cell that includes four Bi and 
four Ni sites, as shown in Fig. [T] d). By investigating several 
sets of parameters and varying them, we find that U B and A 
are most relevant to the valence transition; hence, below we 
show the results by varying these two parameters while fixing 
the others at f N = f B = 1, f BN = 0.5, Un = 3, V B = 0.65, 
Vbn = 1- We have confirmed that the detailed changes of the 
parameters do not alter the following results qualitatively. 

First, we present the ground state properties. Figure [2ja) 
shows the obtained phase diagram on the plane of A and 
U B . There are three main phases with different valence states 
[Figs. [Ajbj-Od)], while two sub-phases appear between them 
[Figs. He) and [2jf ) |. In the uniform-I (II) phase, a major por¬ 
tion of the electrons reside in the Bi (Ni) sites with spatially 
uniform distribution, and no spin polarization is seen in either 
sites as shown in Fig. [2jb) (2d)]. On the other hand, in the 
CO+AFM phase, the average charge density becomes almost 
the same for Bi and Ni, i.e., charge transfer occurs between 
them compared to the uniform phases. In addition, a bipo- 
laronic charge disproportionation occurs in a staggered way 
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FIG. 3: (Color online) A dependences of the charge densities (upper 
panels) and spin moments (lower panels) at each Bi and Ni site for 
/ii - i n - 1, £bn - 0.5, Un = 3, and Vb = 0.65; (a) Ub = —2, 
(b) Ub = 0, and (c) Ub = 2. 

between the Bi sites, as shown in Fig. [2c). The spin moments 
are canceled at each Bi site, whereas those at the Ni sites are 
antiferromagnetic ally ordered. 

To elucidate variations of the electronic states in more de¬ 
tail, we show the A dependences of the charge densities 
( n f^) = (n^ N ^) + (n^ N ' > ) and spin moments (sf N ) = 
(( n it) — The spin moments in the Bi sites remain 

zero in the parameter range we calculated. Figure [2 a) shows 
the results in Ub = —2, where the valence skipping nature of 
Bi is relatively strong. In the large-A region, where the energy 
level of the Ni orbital is substantially higher than that of Bi, 
the electrons occupy dominantly and uniformly the Bi sites, 
and the Ni orbitals are nearly empty; no spin polarization ap¬ 
pears. This electronic state corresponds to the uniform-I phase 
in Fig. [2jb). When A is decreased, the charge densities and 
the Ni spin moments simultaneously jump at A ~ 1.25, with 
showing the spontaneous symmetry breaking in both charge 
and spin channels: the system enters the CO+AFM phase. At 
the transition, the average charge density in the Bi (Ni) sites 
suddenly decreases (increases) by about 0.7. Namely, the va¬ 
lence transition is caused by the charge transfer between Bi 
and Ni. 

As seen in the phase diagram in Fig. [2 a), when Ub is in¬ 
creased, the CO+AFM phase shrinks and shifts to larger A 
region. This is seen in Figs. [2b) and [2c), showing the re¬ 
sults for Un = 0 and Ub = 2, respectively. When Ub = 0, 
the uniform-II phase is realized in the small-A region, where 
the charge density in the Ni sites becomes larger than that in 
the Bi sites, on the contrary to the uniform-I phase. The tran¬ 
sition between the CO+AFM and uniform-II phases is also 
a valence transition. In addition, a narrow intermediate state 
termed sub-II [Fig. [2jf)] appears, in which the average charge 
densities are similar to the uniform-II phase but with a weak 
checkerboard-type CO at the Bi sites. For Ub = 2, the region 



FIG. 4: (Color online) Global phase diagram in the A-Ub-T space. 
Solid and broken lines denote the second and first order transitions, 
respectively. Thin dotted lines are guides for the eye denoting con¬ 
nections between the CO domes. 

of A for the CO+AFM phase is further narrowed and shifted, 
and another intermediate phase termed sub-I [Fig. |2 e )] ap¬ 
pears between the CO+AFM and uniform-I phase, where a 
weak AFM order occurs at the Ni sites. We note that a larger 
Ub results in a smaller amplitude of the valence changes, as 
shown in Figs.[2a)-[2c). 

Next, we discuss the finite-T properties. Figure[4]shows the 
global phase diagram with the T-axis attached to the ground- 
state phase diagram in Fig. |2a). When T is raised in the 
CO+AFM phase, first the AFM order is destroyed, while the 
CO remains; at higher T the CO melts and the uniform phase 
is realized. The uniform-I and -II phases in the ground state 
are continuously connected in the high T region, as they pos¬ 
sess the same symmetry. Consequently, the CO phase shows 
a dome-like structure. The width of the CO dome becomes 
wider for lower T, whose rate is slightly increased below the 
AFM transition temperature. 

The T dependences of the charge densities at each site are 
shown in Figs. [2a) (2c) for several values of A at Ub = —2. 
The CO transition is discontinuous for the large and small val¬ 
ues of A, while it becomes continuous for intermediate A. 
The discontinuous CO transition is accompanied by the large 
amount of the charge transfer between the Bi and Ni sites. 
Below the CO transition temperature, the average charge den¬ 
sities in the Bi and Ni sites become nearly 1. The transition 
between CO and CO+AFM is always continuous, whose criti¬ 
cal temperature is insensitive to A. Figures[2d)-Bf) show the 
results for Ub = 2; all the T scales become smaller compared 
to the case of Ub = —2, and accordingly, the amplitude of the 
CO also becomes smaller. The average charge densities in the 
Bi and Ni sites change even below the AFM transition tem¬ 
perature, which moderately depends on A here, in contrast to 
the case of Ub = —2. 

Through the above analyses, we have found that the valence 
transition occurs between the uniform and CO phases. Here, 
we discuss its origin. The Bi-site electrons have the instability 
toward the bipolaronic CO due to the combination of negative 
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FIG. 5: (Color online) T dependences of the charge densities at each 
Bi and Ni site for f B = fN = 1, Ibn = 0.5, ?7n = 3, and Vb = 
0.65; (a)-(c) 17b = —2 and (d)-(f) Ub = 2 for different values of A 
(indicated in the figures). Note the difference in the T range for the 
two choices of Ub ■ 


Ub and positive Vb■ The instability is most enhanced when 
the Bi band is half-filled; J2i{nf)/N = 1 (the sum is taken 
over the Bi sites and N is number of the Bi sites). This con¬ 
dition is satisfied near the center of the CO dome where the 
CO temperature is maximized, for example in the situation in 
Fig. Ob). On the other hand, apart from the dome center, the 
Bi band filling becomes incommensurate in the high-T uni¬ 
form phase, as shown in Figs. 0 a) and0c). In this case, the 
CO transition is accompanied by the charge transfer between 
the Bi and Ni sites for restoring the commensurability of the 
Bi band. This commensurate locking is the primary cause of 
the valence transition. The AFM instability in the Ni sites due 
to the positive 17 n also contributes to this locking effect, as 
clearly seen in Fig.0d). These considerations lead us to con¬ 
clude that the valence transition is attributed to the commen¬ 
surate locking of the electron filling in both Bi and Ni bands 
driven by the strong electron correlations. 

Finally, let us compare the present results with the experi¬ 
ments in BiNiO.i and related materials. The first-order valence 
transition with the charge transfer from Ni to Bi is observed 
by increasing T and pressure in BiNiO.j. In the present cal¬ 
culation, the discontinuous valence transition driven by T is 
indeed seen, e.g., in Fig. 0 a). The pressure-driven transition 
can also be explained by the transition with increasing the rel¬ 
ative orbital energy A, e.g., in Fig. [3] a), when we assume that 
the pressure affects A. This is reasonable as the volume re¬ 
duction under pressure is mainly due to the contraction of the 
Ni-O perovskite framework in BiNiC >3 11511 , which may in¬ 
crease A via the increase of the Ni 3d energy levels. 

On the other hand, A can be more directly controlled by 
the substitution of the TM cations. Such effects are ac¬ 
tually investigated in B 1 MO 3 (M = V, Cr, Mn, Fe, Co, 


and Ni) raSj and PbMOg (M = Ti, V, Cr, Fe, and 
Ni) [24, j25j], where Pb is another valence skipper favoring the 
oxidation states Pb 2+ (6s 2 ) and Pb 4+ (6s°). In the Bi com¬ 
pounds, the valence state is Bi 3+ M 3+ for M = V - Co, 
while it is Bi 3+ 0 . 5 Bi 5 +o. 5 M 2+ for M = Ni (at ambient pres¬ 
sure) as we have been discussing. Meanwhile, in PbMOs, 
successive variations of such valence state as Pb 2+ M 4+ -H- 
Pb 2+ o. 5 Pb 4 +o.5M 3+ o Pb 4+ M 2+ are observed between V 
and Cr, and Fe and Ni, respectively, where the intermediate 
state Pb 2+ o.5Pb 4+ 0 .5M 3+ involves the CO of Pb 2+ and Pb 4+ 
resembling Bi 3+ o. 5 Bi 5+ o. 5 -ftT 2+ . These experimental results 
suggest that the electrons tend to be accumulated in the A(B) 
cations in the systems with low (high) 3d level, while they are 
uniformly distributed in A and B cations for the intermediate 
case. This tendency and the phase transitions between the dif¬ 
ferent valence states are well understood by the A dependence 
in our phase diagram in Fig. [4j 


In summary, we have presented and analyzed a microscopic 
model for the electronic states in the perovskite oxides includ¬ 
ing the valence skipper as the ,4-site cation, especially focus¬ 
ing on the valence transition in BiNiO j. We have found that 
the valence transition is attributed to the commensurate lock¬ 
ing of the electron filling in the Bi and Ni bands due to the 
electron correlations, and it is sensitive to the relative energy 
of the Bi 6s and Ni 3d levels and the intra Bi-site interac¬ 
tion. Our work provides a fundamental understanding of the 
electronic properties in a series of new perovskite materials 
including valence skippers as the ,4-site cation which have as 
yet been scarcely dealt with. As a mechanism of negative ther¬ 
mal expansion, the ion radius change due to the valence tran¬ 
sition is markedly distinct from those previously discussed, 
which are attributed to lattice vibrations and magnetovolume 
effects lf26ll . We expect that our result will contribute to the 
development of the “next generation” negative thermal expan¬ 
sion materials. 
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